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We investigate the possibility of projecting low dimensional chaos from spatiotemporal dynamics 
of a model for a kind of plastic instability observed under constant strain rate deformation condi- 
tions. We first discuss the relationship between the spatiotemporal patterns of the model reflected 
in the nature of dislocation bands and the nature of stress serrations. We show that at low applied 
strain rates, there is a one-to-one correspondence with the randomly nucleated isolated bursts of 
mobile dislocation density and the stress drops. We then show that the model equations are spa- 
tiotemporally chaotic by demonstrating the number of positive Lyapunov exponents and Lyapunov 
dimension scale with the system size at low and high strain rates. Using a modified algorithm for 
calculating correlation dimension density, we show that the stress-strain signals at low applied strain 
rates corresponding to spatially uncorrelated dislocation bands exhibit features of low dimensional 
chaos. This is made quantitative by demonstrating that the model equations can be approximately 
reduced to space independent model equations for the average dislocation densities, which is known 
to be low-dimensionally chaotic. However, the scaling regime for the correlation dimension shrinks 
with increasing applied strain rate due to increasing propensity for propagation of the dislocation 
bands. 
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I. INTRODUCTION 

It is well known that a number of irregular scalar exper- 
imental signals have features of low dimensional (low-d) 
chaos. These signals often correspond to a spatial av- 
erage over internal degrees of freedom of a spatially ex- 
tended system. In such cases, while the internal degrees 
of freedom are not accessible to experiments, they are 
suspected to be spatio-temporally chaotic, a view sup- 
ported by models that capture the basic features of the 
phenomenon. A good example is the plastic deforma- 
tion of metallic alloys subjected to a constant strain rate 
test. Here, only stress can be measured by the load cell 
placed at one end of the sample and is the spatial av- 
erage of dislocation activity in the sample. Under spe- 
cific conditions, the measured stress-time series, cr(t), ex- 
hibits irregular serrations [1, 2]. Such stress signals have 
been shown to be low-d chaotic [3-5]. A second exam- 
ple is the crackling audible noise, called acoustic emission 
(AE) commonly experienced during peeling of an adhe- 
sive tape. The AE signals have also been shown to be 
low-d chaotic [6] , which is controlled by the rugged nature 
of the peel front. These two examples suggest that spa- 
tial averaging some how projects low-dimensional chaos 
out of spatio-temporal dynamics. Several candidates for 
low-d chaos could be cited, including the irregular volt- 
age fluctuations measured across a sample of a charge 
density wave compound (CDW), attributed to pinning 
and unpinning of CDW [7], and the fluctuations in the 
position and intensity of a light beam passing through a 
turbulent medium [8]. 



* ritupan@mrc.iisc.ernet.in 
t garani@mrc.iisc.ernet.in 



The purpose of this paper is to examine the possi- 
bility of projecting low-d chaos from a spatiotemporal 
chaotic system and if so, under what conditions? This 
question will be addressed within the context of a model 
for a plastic instability, called the Portevin-Le Chatelier 
(PLC) effect [9-12]. The PLC effect or jerky flow refers 
to the irregular stress-strain curves associated with the 
observed heterogeneous deformation when specimens of 
dilute metallic alloys are deformed in a window of strain 
rates and temperatures. While dislocation activity in the 
sample is not accessible to experiments, dislocation bands 
seen on the surface of the sample [2, 13], measurement 
of acoustic emission studies [14], and studies on micro- 
crystals [15], strongly suggest that dislocation dynamics 
is intermittent. However, relating the irregular nature 
of the scalar signal to the dislocation dynamics has re- 
mained a difficult task due to lack of dislocation based 
models even though there is a general consensus that dis- 
location dynamics is intermittent, a view supported by 
dislocation dynamics simulations as well [16, 17]. 

The basic mechanism attributed to the stress serra- 
tions in the PLC effect is the collective pinning and un- 
pinning of dislocations [1, 2]. Each stress drop is as- 
sociated with the formation and often the propagation 
of dislocation bands. Three distinct types of dislocation 
bands, namely, the static uncorrelated type C, hopping 
type B and continuously propagating type A are seen 
with increasing strain rate. The average stress drop size 
is large for the type C band, decreasing with increasing 
strain rate as we encounter the type B and A bands [2]. 

Several of these generic features of the PLC effect are 
captured by the Ananthakrishna (AK) model [2, 9-12]. 
One prediction specific to the model relevant for the cur- 
rent study is the (low-d) chaotic nature of the stress drops 
for low strain rates [18]. The prediction has been subse- 
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quently verified using experimental signals from single 
and polycrystals [3-5]. Thus, the model provides a plat- 
form to verify the conjuncture that the stress signals are 
low-d chaotic while the model equations are spatiotempo- 
rally chaotic. Using the AK model, we show that at low 
strain rates, the irregular stress-time series can be un- 
ambiguously identified as low-d chaotic while the model 
equations are spatiotemporally chaotic. We show that 
the low-d chaotic nature is entirely due to the one-to- 
one correspondence between the isolated bursts of mobile 
density and the stress drops in this domain of applied 
strain rates. The low-d chaotic feature breaks down for 
higher applied strain rates which can be directly corre- 
lated with the propensity for the dislocation bands to 
propagate. 



II. THE ANANTHAKRISHNA MODEL FOR 
SERRATED FLOW 

The basic idea of the AK model is that most generic 
features of the PLC effect such as the existence of the in- 
stability in a window of strain rates, the negative strain 
rate sensitivity of the flow stress, and the three types of 
bands (C, B, and A) emerge from the nonlinear interac- 
tion of a few dislocation populations assumed to repre- 
sent the collective degrees of freedom of the system [9- 
12]. The model uses three types of dislocation densities: 
the fast mobile p m (x,t), the immobile pi m (x,t), and the 
Cottrell type p c (x,t) corresponding to dislocations dec- 
orated by solute atoms [9]. (The standard AK model 
equations take into account the hardening process that 
is relevant for describing the physical phenomenon [9-12]. 
However, since we are concerned here with the spatiotem- 
poral dynamics in a stationary situation, wc ignore the 
hardening term that also helps to attain the stationary 
state faster.) The scaled equations for the three densities 
and the scaled stress <f> are given by 



dt 
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The term i>oPm m Eq. (1), refers to the formation of 
dipoles and other dislocations locks, p m Pim refers to the 
annihilation of a mobile dislocation with an immobile 
one, and pi m represents the reactivation of the immobile 
dislocation due to stress or thermal activation. ap m rep- 
resents the immobilization of mobile dislocations due to 
aggregation of solute atoms to dislocations. Once a mo- 
bile dislocation starts acquiring solute atoms we regard it 



as Cottrell- type dislocation p c . As more and more solute 
atoms aggregate, they eventually stop, and are consid- 
ered as immobile dislocations pi m . This term then acts 
as a source term in Eq. (3). <j) m p m represents the rate 
of multiplication of dislocations due to cross slip. This 
depends on the velocity of the mobile dislocations taken 
to be V m ((j>) = <p m , where m is the velocity exponent. 
Within the scope of the model, cross slip that allows 
dislocations to spread into neighboring spatial locations 
gives rise to diffusive coupling [the last term in Eq. (1)]. 
These equations are coupled to Eq. (4) , which represents 
the constant strain rate deformation experiment. In Eq. 
(4), e a is the scaled applied strain rate, d the scaled ef- 
fective elastic modulus of the machine and the sample, 
and I the dimensionless length of the sample. The scaled 
constants, a, bo, and c refer, respectively, to the concen- 
tration of solute atoms slowing down the mobile disloca- 
tion, the thermal and athermal reactivation of immobile 
dislocations, and the diffusion rate of solute atoms to the 
mobile dislocations. The drive parameter is the applied 
strain rate e a with respect to which the different types 
of bands and the associated serrations are observed. The 
instability range is found in the interval 20 < i a < 1625 
for a system size N — 100. 

Equations (1—4) are discretized on a grid of N points 
and solved using an adaptive step-size differential equa- 
tion solver (MATLAB "ode23"). The initial values for 
the densities are taken to be uniformly distributed along 
the sample. However, as the long-term evolution does 
not depend on the initial values, steady-state values have 
been used. As dislocation bands can not propagate into 
the grips, the boundary values Pi m {j, t); j = 1 and N are 
taken to be two orders higher than p% m {j),j — 2, ..N — 1. 
Further, we impose p m (j,t) — p c (j,t) — for j — 1 and 
N. The results reported are for a = 0.8, b = 5 x 10 -4 , c = 



0.08, d = 6 x 10" 



3,D = 0.25. The system size 



N = 20 — 150 is used depending on the property ad- 
dressed. After discarding the initial transients, 4>(t) is 
sampled at time intervals 5t = 0.5. 



III. SPATIO-TEMPORAL PATTERNS AND 
RELATION TO STRESS SERRATIONS 

The results of the model relevant for the present pur- 
pose are the spatiotemporal patterns reflected in the na- 
ture of the dislocation bands and their connection to the 
changes in the nature of stress serrations, observed with 
increasing applied strain rate. To understand the origin 
of the instability and the nature of the spatial patterns, 
let us first examine Eqs. (1) and (4) more carefully. We 
first note that <f> = o j a y , where a is the unsealed stress at 
time t and a y is the yield stress. From Eq. (1), it is clear 
that p m increases abruptly only when (j> exceeds unity 
(corresponding to the unsealed stress exceeding the yield 
stress) or equivalcntly, there is a threshold for nucleation 
of an isolated burst of the mobile density p m . However, 
a stress drop can only occur when the space-averaged 
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plastic strain rate e p — - ^ f Q p m (x,t)dx exceeds the 
applied strain rate e a . 

As mentioned earlier the model reproduces all the three 
types of bands (C, B, and A) observed in experiments 
with increasing e a [2, 12]. For a range of low e a val- 
ues from 20 to 50, randomly nucleated dislocation bands 
are seen in the form of isolated bursts of p m . From the 
above discussion, we know that an isolated burst of p m 
can only be nucleated if the scaled stress c/> increases be- 
yond unity. A typical space-time plot of the randomly 
nucleated bands is shown in Fig. 1(a) for i a = 40. The 
corresponding stress-time curve is nearly regular as illus- 
trated in Fig. 1(b). Furthermore, in this range of e a , 
usually a single burst of p m contributes to the total plas- 
tic strain rate e p — - ^ J p m (x,t)dx as the magnitude 
of p m (x, t) at other spatial locations is insignificant. This 
implies that whenever the space-averaged plastic strain 
rate overshoots the applied strain rate e a , a stress drop 
occurs [see Eq. (4)]. Thus, there is a one-to-one corre- 
spondence between the burst of mobile density p m and 
the stress drop. However, as we increase e a , a new burst 
of p m is formed even before the previous burst dies off giv- 
ing the impression of hopping character (type-B band). 
This is clear from the plot of the hopping-type partial 
propagating band shown in Fig. 2(a) for e a = 90. This 
also means that the stress required to nucleate a fresh 
burst of p m ahead of the previous one even before it dies 
off, is less than that required to nucleate an isolated burst 
of p m . This in turn implies that the amplitude of the 
serrations would be smaller in regions where propagat- 
ing bands are seen. This relationship between the prop- 
agative nature of the bands and small-amplitude stress 
drops is illustrated in Figs. 2(a) and 2(b). The figure 
identifies the propagative regime with the corresponding 
stretch of small amplitude stress serrations [marked by a 
corresponding set of arrows in the Fig. 2(a) and 2(b)]. 
As we increase e a , the extent of propagation increases, 
concomitantly, longer stretches of small-amplitude serra- 
tions interrupt the otherwise large-amplitude serrations. 
The spatial correlation increases with e a until the bands 
propagate fully with numerous small stress drops. The 
corresponding stress-time plot would have mostly small- 
amplitude stress drops with very few large stress drops. 
A typical plot of a fully propagating band along with the 
associated stress-time plot is shown in Figs. 3(a) and 
3(b). (see Refs. [2, 12]). 
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FIG. 1. (Color online) (a) Randomly nucleated type-C bands 
for e a — 40. (b) The corresponding stress-time plot. 
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FIG. 2. (Color online) (a) Hopping type-B bands for e a = 
90. (b) The corresponding stress-time plot. For the sake of 
clarity, a short segment of the space-time plot is displayed. 
The correspondence between the propagating nature of the 
band marked by the arrows in (a) and the stress-time series 
in (b) is displayed. 



IV. SPATIO-TEMPORAL DYNAMICS OF THE 
MODEL 

The above spatio-temporal dynamics has been quan- 
tified by calculating the Lyapunov spectrum from the 
model equations (using the Bennettin's algorithm [19]) 
for a range of values of e a where the randomly nucleated 
bands and partial propagative bands are seen. (First 
10, 000 points have been ignored and the spectrum is cal- 
culated using the next 10, 000 time steps.) We have stud- 



ied the system size dependence of the number of positive 
exponents and the Lyapunov dimension D^. Both 
scale linearly with N for low strain rates where randomly 
nucleating bands are seen and also at high strain rates 
where fully propagating bands are seen. For illustration, 
plots of Dl verses N is shown in Figs. 4(a) and 4(b), 
for e a — 40 and e a — 240, respectively. However, for a 
range of values of e a beyond 60 where partial propaga- 
tive bands are seen [that also coincides with the region 
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where we do not find converged values of the correlation 
dimension, £>2s (f, d) , as we shall soon show] , we find two 
distinct slopes in the Dl — N plot, one for small values 
of N and another for large values of N. This feature can 
be seen even for e a — 60 where partial propagation be- 
gins. A typical plot of Dl — N is shown in Fig. 5(a) for 
e a = 120. (The larger slope for small N values is ~ 0.67 
and the smaller slope for large N is ~ 0.287.) As we 
increase the strain rate, the range of values of N corre- 
sponding to the larger slope decreases; eventually a single 
slope is seen for e a — 240 [Fig. 4b]. A similar behavior 
is also seen for rit. 

To understand the underlying reasons for the two slope 
nature of Dl versus N, we have examined the nature 
of the spatiotemporal patterns when the system size is 
increased from small values of N to large values. Plots 
of spatiotemporal patterns for e a = 120 for N = 100 and 
N = 20, are shown in Fig. 5(b) and 5(c), respectively. As 
can be seen from the two plots, while the bands propagate 
fair distances for large system size N = 100, for small 
system size N — 20, one finds well-separated bursts of 
p m (in the time domain) at any given time, although 
there is a visual impression of propagation. However, the 
direction of apparent propagation [in Fig. 5c] is opposite 
to that for N = 100 [Fig. 5b]. Indeed, we find that the 
dynamics is altered from partial propagative nature for 
large system size (say for N = 100) to one of burst type 
when the system size is smaller than N = 50. (Note that 
the burst-type pattern is seen for e a < 60 when N is 
large.) Similar altered dynamics for small system sizes is 
seen for 60 < e a < 200. 



FIG. 4. (Color online) Lyapunov dimension Dl as a function 
of TV (a) for e a = 40 (randomly nucleating bands) with slope 
w 0.35 and (b) for e a = 240 (continuously propagating bands) 
with slope ss 0.37. 



V. TIME SERIES ANALYSIS OF STRESS 
SIGNALS 

Now consider analyzing the stress-time series that is re- 
lated to the spatial average of the dislocation activity in 
the sample [Eq. 4] . Consider a stress-time series of length 
M in units of St defined by {<t>(k),k = 1,2, ••• ,M}. 
Then, the reconstructed attractor is defined by the cIe 

dimensional vectors ^- = {<j>(k), (f>{k+r), ■ ■ ■ ,(j)[k + (d,E — 
1)t]}; k = 1,2, • • • , [M - (d E — l)r], where r is the delay 
time. The chaotic nature of the attractor is quantified by 
establishing the existence of finite correlation dimension 
and a positive Lyapunov exponent. 

The most popular method for calculating the correla- 
tion dimension is the Grassberger-Procaccia (GP) algo- 
rithm [20]. The method calculates the correlation inte- 
gral defined as the fraction of the pairs of points £j and 
£j whose distance is less than a specified value r, that is, 
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where 0(...) is the Heaviside step function and M p the 
number of vector pairs used [20]. A window is imposed 
to exclude temporarily correlated points. The correlation 
dimension is defined as 



D 2 {r,d E ) 
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FIG. 5. (Color online) (a) Lyapunov dimension Dl as a func- 
tion of N for e a = 120. Note the two distinct slopes. One 
for larger system sizes (N > 50) with slope ~ 0.287 (dashed 
line) and the other for smaller system sizes with slope ~ 0.67 
(continuous line) (b) Space-time plot of the mobile dislocation 
density for the same e a and N = 100. (c) Space-time plot of 
the mobile dislocation density for e a = 120 and N = 20. 



This limit is seldom reached. Instead, one usually finds 
a reasonably large scaling regime in In r at intermediate 
scales, where the slope D 2 (r,d,E) converges to a finite 
value D 2 , which is taken to be the correlation dimension 
of the attractor. The method has been successfully ap- 
plied for analyzing several (mostly) low-d attractors [21]. 
However, the validity of the GP algorithm for high di- 
mensional attractors has been questioned since exponen- 
tially longer time series are required [22-24]. Generally, 
two types of errors limit the confidence in the estimates of 
D 2 , namely the poor statistics at small length scales and 
the underestimation of D 2 (r,dE) at length scales com- 
parable to the attractor size. The latter arises when the 



FIG. 6. (Color online) Plot of D 2 s(d E ,r) as 

a function of Inr for the sum of two independent 
x components of the Lorenz model for (c£e,t) = 
(5, 40), (6, 32), (7, 27), (8, 23), (9, 20), (10, 18). Normalization 
used is D 2 (r, 2)/2 for r = 250. 



reference point is close to the edge of the attractor as 
zero contribution arises for length scales r beyond the 
attractor's edge. This systematic error is shown to be 
important in characterizing high dimensional attractors 
including those of spatially extended systems [25, 26]. 
Bauer et al. [25] suggested a method of compensating the 
contribution arising from finite size of the attractor by 
normalizing the local slope D%(r,d E ) by the slope of an 
equivalent random attractor A2r(?">^£) = c ""^^' dE ^ ■ 
Thus, the "true" correlation dimension is redefined as, 



D 2s (r,d E ) = 



D 2 {r, di 



dlnC2(r, dE) 



D2r(r,d E ) d E dlnC 2r (r,l) 



(7) 



where C 2r (r,dE) — C 2r (r,l) dE has been used [25]. Here 
we use C 2 (r, 2) 1 / 2 instead of C 2r (r, 1) and define 

n , , s gafc d E ) ^ dlnC 2 {r, d E ) , a . 
°^ dE) = = 2 dlnC 2 (r,2) ' (8) 

Using C 2 (r, 2) 1 / 2 not only serves to correct for the finite 
size of the attractor as does C 2r (r, 1), it also includes the 
contribution from finite delay time. The converged value 
of D 2s (r,dE) over a fair range of Inr is taken to be D 2 . 
We further use constant window length (d E — 1)t that 
maximizes the scaling regime [27]. The Lyapunov spec- 
trum for the time series is computed using the standard 
Eckmann's algorithm [28]. 

We first check the applicability of the algorithm to 
high dimensional attractors with dimensions, say four 
or five. For this, we consider a time series obtained 
by summing two independent x-component time series 
of the Lorenz model. The first data set is for parame- 
ter values a — 10, r = 28,6 = 8/3 and the second is for 
a = 16, r = 40,6 = 8/3. The total number of points 
used is 10 5 . We have kept (d E — 1)t = 160, an op- 
timum value of the window length that maximizes the 
scaling regime to calculate the C 2 (r,dE) for (dE>T~) = 
(5, 40), (6, 32), (7, 27), (8, 23), (9, 20), (10, 18). The slope 
D 2 (r, ds) has been normalized with respect to D 2 (r, 2) /2 
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FIG. 7. (Color online) Plot of D2s(r, <fe) as a function of Inr 
for the stress-strain series for i a = 25 for 1.1 x 10 s points. The 
values of (ds,r) are (4, 16), (5, 12), (7, 8), (9, 6). Normalization 
used is D 2 (r, 2)/2 for r = 30. 



for t = 65. Figure 6 shows the plot of D2 S (r, ds) for this 
case. It is clear that there is nearly five orders of scal- 
ing regime, which is comparable to the scaling regime 
obtained using an adoptive box counting method where 
10 7 points have been used [29] . The value of Di obtained 
is ~ 3.8. Note that 2L>2 = 4.12 is only the upper limit. 
The above algorithm works very well for low-d attractors 
such as the Lorenz model (where we find six orders of 
scaling regime with just 10,000 points). 

Having demonstrated the spatiotemporally chaotic na- 
ture of the model equations, we now examine if the stress- 
time series has the required features of low-dimensional 
chaos using the above algorithm. Here we note that the 
stress rate 4> depends linearly on the plastic strain rate 
which in turn involves the spatial average of mobile dis- 
location density. Thus, the scalar signal <p, though a 
dynamical variable, which determines the rate of multi- 
plication of mobile density [first term in Eq. (1)], is an 
appropriate quantity for carrying out time series analysis. 
We have calculated D2 from the stress-time series <p(t) in 
the range 20 < e a < 50 corresponding to the randomly 
nucleating bands using 1.1 x 10 5 points for N = 100. 
A plot of D2 S (r,d,E) as a function of Inr for i a = 25 is 
shown in Fig. 7, keeping (d,E — l)r = 48 (cIe = 4, 5, 7 and 
9) [27]. At least five orders scaling regime in Inr is clear 
and D2 = 2.15 ± 0.03. D2 increases with e a marginally 
until 50. Beyond this value, we find the scaling regime 
shrinks. Indeed, even for e a = 60, Z^s^^e) increases 
steadily but slowly, for small length scales as is clear from 
Fig. 8(a). Even when the propagation distance is not too 
large we find practically no scaling regime as illustrated 
in Fig. 8(b) for e a = 90. The nonconstancy of D2 S (r, <Ie) 
for e Q > 50 is precisely the region of the partially propa- 
gating dislocation bands (see Fig. 2(a)). 

Using the Eckmann's algorithm, we have calculated 
the Lyapunov spectrum for the time series in the interval 
20 < e a < 50. Fig. 9 shows the spectrum for e a = 25 
for d,E — 6. ( Note that d,E — 6 corresponds to 2D 2 + 1. 
However, Ding et al. [30] have shown that it is adequate 
to use embedding dimension cLe > D2 rather than using 
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FIG. 8. (Color online) (a) Plot of D2 S versus In r for the stress- 
time series for e a = 60. The values of (<1e,t) are such that 
(cLe — 1)t = 27. Normalization used is D2(r,2)/2 for r = 
18. (b) Plot of D2S versus Inr for the stress-time series for 
e a — 90. The values of (ds, r) are such that (ds — l)r = 24. 
Normalization used is D2(r, 2)/2 for r = 16. 



a (Le larger than 2D 2 + l as originally suggested.) The ex- 
istence of a good zero exponent and a positive exponent 
is clear and the Lyapunov dimension Dl = 2.78. Similar 
results are obtained for 20 < e a < 50. D L ranges from 
2.6 — 2.8. The minimum degrees of freedom required for 
a dynamical description is then four, which is also the 
dimension of the bare (space independent) AK model. 
Here, we note that as we increase the embedding dimen- 
sion d,E from D2 + 1 to 2Z?2 + 1, the value of the positive 
exponent remains nearly constant at 0.075, the zero 
exponent remains close to zero (~ 10 -3 ) and the first 
negative exponent also remains nearly constant. These 
three exponents are nonspurious. For e a > 60, it is not 
meaningful to calculate the Lyapunov spectrum due to 
lack of convergence of -D 2s (r, tig). 

We have also carried out surrogate (phase-randomized) 
data analysis of the stress-time series for e a = 25 — 50. 
We find that D2 S (r,dE) increases with the embedding 
dimension. This provides an additional support that the 
time scries has features of low dimensional chaos. 



VI. SUMMARY AND CONCLUSIONS 

In summary, we have demonstrated that low- 
dimensional chaos is detected in the stress signals for the 
low-strain-rate regime even as the model equations are 
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FIG. 9. (Color online) The Lyapunov spectrum for cLe = 
6 and r = 2 for the stress-strain series for t a = 25. The 
Lyapunov dimension is 2.78. 



spatiotemporally chaotic. The scaling regime for the cor- 
relation dimension shrinks beyond e a — 50 where the ten- 
dency for propagation increases with the applied strain 
rate. The number of positive Lyapunov exponents and 
the Lyapunov dimension for the model equations scale 
with the system size for the low strain regime of ran- 
domly nucleated bands. These two results appear to be 
mutually inconsistent if one goes by the general belief 
that a scalar time series obtained from a spatially ex- 
tended system should in principle contain information 
about the full attractor of the system. However, the im- 
pressive (five orders in In r at low applied strain rates) 
scaling regime for D 2 for the stress-time series strongly 
supports the low-dimensional chaotic nature of the stress 
signals. Interestingly, the value of D 2 for the stress-time 
series obtained from the bare (space-independent) AK 
model equations is also ~ 2.25. Furthermore, the Lya- 
punov spectrum calculated from the embedded time se- 
ries shows a positive and a good zero exponent, which 
confirms the low-d chaotic nature of the stress signal. 

On the other hand, the scaling regime for the correla- 
tion dimension shrinks for higher strain rates (e > 50) 
once the bands acquire the propensity to propagate. In 
the partially propagative regime the Lyapunov dimen- 
sion D L exhibits two slopes as a function of the system 
size N, one for the small sizes (N < 50) and another for 
larger sizes while we find Dl again scales with N for the 
fully propagating bands at high strain rates as for the low 
strain rates. We have shown that the two-slopes nature 
is due to the altered spatiotemporal patterns when the 
system size is reduced from those manifesting for large 
system sizes. 

Having demonstrated that the stress signals have all 
features of low-dimensional chaos for the region of strain 
rates 20 < e a < 50, a natural question is as follows: 
How is the low-d chaos projected from spatiotemporal 
chaotic system? The answer lies in the nature of spa- 
tiotemporal patterns seen at low strain rates. Plots of 
p c {x,t) and pi m (x,t) corresponding to p m (x,t) shown 
in Fig. 1(a) for e a = 40 are displayed in Figs. 10(a) 
and 10(b), respectively. Since the relative magnitudes of 



the three densities are substantially different, this fea- 
ture is better captured by a snap shot representing the 
magnitudes of the three densities (at an arbitrary time) 
as shown in Fig. 10(c). As demonstrated earlier, there 
is a one-to-one correspondence between the stress drops 
and bursts of p m . Furthermore, on comparing p m (x,t) 
shown Fig. 1(a) with p c (x,t) [Fig. 10(a)], it is clear that 
these two densities are in phase and are localized to the 
same spatial extent of few sites with large peak heights 
of p m (x,t) ~ 1000 and p c (x,t) ~ 500 respectively. In 
contrast, the range of values of p im (x,t) is two to three 
orders lower than p m and p c as can be seen from Fig. 
10(c). Moreover, the range of pi m is small with values 
between 0.5 — 3. [Note that larger values of pi m are for 
sites where p m (or p c ) is small and vice versa, as can be 
seen from Fig. 10(c).] Moreover, we note that the spa- 
tial dependence in the model comes entirely from Eq. (1). 
Using these features, we shall now show that the spatial 
averaging process does project low-dimensional chaos. 

To do this, we first note that since the contribution to 
the plastic strain rate [i p = - ^ J Q p m (x,t)dx] comes 
from the single burst, it is natural to use the space- 
averaged dislocation densities to illustrate the projection 
process. Second, since both p m and p c are localized to a 
few sites with large peak heights, they can be represented 
by Gaussian-like functions. Now consider integrating Eq. 
(!)• 



dp„ 
dt 



dx = 



J dx[- b p 2 m - PmPim + Pim 

D<p m {t) d 2 Pm , 



apr, 



+<f> m pm + 



dx 2 



(9) 



Then, define J p m (x,t)dx/l = p m (t) with similar defini- 
tions for p c {x, t) (p c (t)) and p lm {x, t) {p im {t)). The linear 
terms in p m (x,t) pose no problems. Using a Gaussian 
representation for p m (x,t), the first term J dxp 2 m {x,£)jl 
in Eq. (9), can be integrated to give J dxp 2 m (x,t)ll = 
-^j^p 2 m {t) where a is the variance of the Gaussian distri- 
bution. 

Consider evaluating J p m (x,t)pi m (x,t)dx/l. This in- 
volves integrating the sharply peaked p m (x,t) with the 
weight factor pi m (x,t). Noting that pi m (x, t) varies in 
a narrow range of 0.5 — 3, and p m (x,t) is localized to 
a few sites, we may approximate J p m (x,t)pi m (x,t)/l sa 
qp m {t)p im {t) where / p im (x,t)dx/l = p tm (t) with a pref- 
actor q to account for the average weight factor arising 
from pi m (x,t) to the integral. This approximation has 
been verified numerically for various intervals of time. 
The factor q is around 2. 

Now consider evaluating the last term in Eq. 9. Noting 
again that the range of pi m (x, t) is nearly constant of the 
order of unity, we can replace it with rpi m (t), where r is 
a scale factor. Then, using a Gaussian centered around 



tions can be integrated to give 
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FIG. 10. (Color online) (a) Space-time plot of p c (x,t) for 
e a — 40 for the same time interval as for p m (x,t) shown in 
Fig. 1(a). (b) The corresponding space-time plot for p im (x, t). 
Note that for p% m {x, t) is out of phase with p m and p c for all 
sites, (c) A snap shot of p m {x,i), p c (x,t) and Pi m (x,t) at an 
arbitrary time. 



some site, we get 

D(b m (t) d 2 p m (x,t) 



I = 



dx 



Pim(x,t) dx 2 

When the p m burst is near the boundary, this term would 
be nonzero. However, such events are rare and for all 
purposes, they may be ignored. Thus we have 



dpn 



dt 



= -b'oPm - QPmPim + Pirn - aPm + 4*™ Pm- (H) 



Here, b' = b ^. 

It is straight forward to show that the other three equa- 



dpi m 

dt 

dp c 

dt 
dcj){t) 

dt 



b' (b'oP m - qPmPim. - Pirn+O-Pc), (12) 



c(Pm ~ Pc), 

d[e a -p m {tW n (t)]. 



(13) 
(14) 



Thus, for low strain rates, Eqs. (1) — (4) reduce to a set of 
coupled ordinary differential equations for the space av- 
eraged densities given by Eqs. (11) — (14), with renormal- 
ized co-efficients. These equations have the same form as 
the bare AK model equations that have been shown to 
be chaotic (see Refs. [2, 18]). Indeed, the value of D 2 
for the stress-time series obtained from the bare model 
turns out to be D 2 — 2.25±0.05 using just 30, 000 points. 
Clearly, the above procedure works as long as there is a 
one-to-one correspondence between the bursts of p m and 
stress drops. However, this correspondence breaks down 
even for a small extent of propagation. Basically, during 
propagation, the successive bursts of p m ahead of the pri- 
mary burst do not have any specific relation to the stress 
drop. 

The fact that the zeroth order Fourier component cap- 
tures the low-dimensional chaotic nature of the stress sig- 
nal is possibly suggestive of separation of time scales. 
Here, we note that the burst time scale of p m (x,t) is 
short and the spatial averaging process maps the burst 
time scale to the time scale of the stress drop. Then, 
the one-to-one correspondence between the burst in p m 
and stress drop leads to the low-d chaotic character of 
the stress reflected in the constancy of D2 S (r, c?_e) for all 
but small length scales. However, in principle, the influ- 
ence of higher Fourier components on the zeroth Fourier 
component must be reflected in some form. The effect of 
other degrees of freedom has been represented as noise to 
explain the generally increasing trend of the correlation 
dimension for small length scales in studies on coupled 
maps [31]. Indeed, the increasing trend of D2 S (r, fig) for 
small scales seen in Fig. 7 is perhaps a reflection of the 
high-dimensional nature of the full system. 

Now consider the Lyapunov spectrum obtained from 
the two methods. The main contribution to the diver- 
gence of orbits of the full set of model equations comes 
from the p m {x,t) and p c (x,t) due to the large range of 
values while both pi m and stress contribute minimally. In 
contrast, the contribution to the divergence of the orbits 
in the embedded space of </>(t) comes mainly from the 
rapid changes occurring during stress drops. Since the 
underlying mechanisms and magnitudes of the changes 
contributing to the Lyapunov spectrum in the two cases 
are very different, the results are not necessarily mutually 
inconsistent. 

We believe that projecting low-dimensional chaos from 
spatiotemporal chaotic dynamics should hold at least in 
situations where there is a one-to-one correspondence be- 
tween the abrupt variation of a scalar time series (which 
is some kind of average over spatial degrees of freedom) 



9 



and localized excitations of the internal degrees of free- 
dom. This could also arise due to synchronization of 
a certain region of spatial elements. Attempts to ver- 
ify this conjecture are underway. Detailed investigations 
arc in progress to understand the relationship between 
the stress signals and the spatiotemporal dynamics of 
the model equations, and quantifying stress signals for 
higher strain rates where partial propagating bands are 



seen. 
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